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Abstract. We argue that the Fermi-Hubbard Hamiltonian describing the physics of ultracold 
atoms on optical lattices in the presence of artificial non-Abelian gauge fields, is exactly 
equivalent to the gauge theory Hamiltonian describing Dirac fermions in the lattice. We show 
that it is possible to couple the Dirac fermions to an "artificial" gravitational field, i.e. to 
consider the Dirac physics in a curved spacetime. We identify the special class of spacetime 
metrics that admit a simple realization in terms of a Fermi-Hubbard model subjected to an 
artificial SU(2) field, corresponding to position dependent hopping matrices. As an example, 
we discuss in more detail the physics of the 2+1D Rindler metric, its possible experimental 
realization and detection. 
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1. Introduction 

The studies of ultracold quantum matter in artificially designed external gauge fields is 
one of the most rapidly developing areas of physics of ultracold atoms Q] [2] [3] |4j |5). 
Originally, these studies arose from investigations related to the response of superfluids, 
such as Bose-Einstein condensates (BEC), to rotation (cf. J6j |7)). On one hand, rotation 
induces quantized vortices and/or vortex Abrikosov lattices (8). On the other, its effects 
are equivalent to those of an artificial constant (Abelian) magnetic field. The latter analogy 
has led to the idea of realizing strongly correlated quantum liquids, such as the celebrated 
Laughlin liquid of the fractional quantum Hall effect (FQHE) (cf. J9J) by means of rapid 
rotation iflOl [TP . Unfortunately, reaching the FQHE regime with rotation is experimentally 
very difficult; it has been achieved recently, but in a system of only 1 < N < 10 atoms 
in rotating microtraps (lattice site potential wells) fPH IT3ll . Several researchers proposed, 
thus, alternative approaches involving for instance laser-induced gauge fields that employ 
dark states in 3-level systems |fl4l . or laser induced gauge fields acting on atoms confined in 
an optical lattice lfl5l [161 [l7ll . Other approaches concerned rotating optical lattices |[T8l[T9l , 
or interactions of lattice atoms embedded in a rotating BEC ll20l . Interestingly, the proposals 
involving lasers can be relatively straightforwardly generalized to particles possessing internal 
"color" states subjected to artificial non- Abelian gauge fields l2"Tl l22l . 

In the last two years there has been a large number of works reexamining these ideas 
and proposing experimental realization within the reach of the current state of the art. The 
NIST group employed an approach similar to lfT4l and used Raman (Bragg) transitions in 
Sodium to realize experimentally the first non-zero constant vector (corresponding to zero 
"artificial" magnetic field) |23l . Note that when the gauge symmetry is provided by Nature 
only gauge invariant observables are physical. In the case of artificial gauge potentials the 
situation is in principle different. Gauge potential are controlled by the experimenters and the 
measurements of quantities that depend on the choice of gauge are possible; in another words, 
the very process of the creation of the gauge potential is not gauge invariant, even though the 
resulting Hamiltonians are (for discussion see 112411251 and l26l ). NIST group was also able 
to generate vortices using the same scheme with a potential configuration corresponding to a 
non zero artificial magnetic field ll27ll . Several practical extensions of the scheme of iTBll were 
discussed in |f28l[29l. New schemes were proposed for alkali and earth-alkali atoms in optical 
lattices employing superlattice techniques l30l 31), and on atomic chips 1321 . Very recently, 
the creation of spin-orbit couplings, a special instance of synthetic non-abelian fields, was 
reported [33]. 

Artificial non Abelian gauge fields are particularly interesting, since they provide a 
natural framework to simulate relativistic physics of the Dirac equation. Artificial Dirac 
physics has been in recent years at the center of interest in condensed matter in the context of 
studies of the amazing properties of graphene [33[35][3S[37][3g[3l|13@^ In 
the case of graphene the Dirac points in the dispersion relation appear due to the geometry of 
the underlying hexagonal lattice. This idea can be carried over to cold atoms |46l ; a hexagonal 
optical lattice (OL) with spinor bosons has been recently realized experimentally 1471 . Also, 
other systems not relying on a lattice have been proposed to emulate Dirac particles l48l . 
Other systems where Dirac physics plays a role include superfluid 3He-A |49ll50l (where the 
analogy to Particle physics can be extended to include gauge interaction and Standard model 
phenomena, for instance see (3T)), trapped ions [52,]53][54][55][56], orFermi-Bose mixtures 
El- 

In the case of artificial non-Abelian laser-induced fields (ANALF) the connection to 
Dirac physics was pointed out in l58l . The Dirac equation is responsible for the anomalous 
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integer quantum Hall effect in artificial SU(2) fields on a 2D square lattice, and topological 
quantum phase transitions on 2D hexagonal lattices, as well as FQHE |59l - These situations 
correspond to relativistic physics in 2+ 1 dimensional spacetime [60. 61 , 62, 63 1. Very recently, 
several proposals were formulated for the creation of ANALF in 3+1 dimensional (3D spatial) 
allowing for simulations of Wilson fermions l64l and axion ll65ll quantum electrodynamics 
with ultracold atoms l66l [311 . 

So far, most of the proposals have dealt with constant non Abelian fields strengths (i.e 
constant Wilson loops), staggered Abelian gauge fields 116711681 . or, in the rare cases, fields 
that are linear in the spacial coordinates 11221 |69l . The Dirac equation, resulting in some 
of these situations, corresponds to a Dirac equation in a flat Minkowski space. The crucial 
ingredient of such studies was based on the analysis of the dispersion relation between energy 
and quasi-momentum, as in graphene (45]. Energy bands touch in isolated singular points, 
called Dirac points. In the vicinity of these points, the dispersion relation linearizes, and a 
cone is formed. Dirac physics can then be realized for fermions adjusting the Fermi surface 
to include the Dirac point, and considering low lying energy excitations. 

In this paper we take a different perspective on this issue, employing standard concepts 
from Hamiltonian lattice gauge theory (HLGT) ITUl . We argue that using artificial non- 
Abelian fields in lattices it is possible to simulate with cold atoms a Dirac spinor in the same 
way as it is done in HLGT. In another words, the Hubbard Hamiltonian describing the physics 
of atoms in artificial non-Abelian fields in lattices is exactly the HLGT version of the Dirac's 
Hamiltonian. That is to say, in the scheme we consider, all excitations of the fermion field in 
the lattice are Dirac-like, not only the low lying ones. The advantage of this point of view is 
that it allows one to couple the simulated Dirac-fields to external fields, or quantum fields in 
a straightforward way, which is not so clear in artificial graphene-like systems. 

In particular, we will show that it is possible to consider coupling Dirac fermions 
to an artificial gravitational field, that is to consider the Dirac equation in curved space. 
We will identify and focus on a special class of space-time metrics that admit a simple 
formulation of the Dirac lattice Hamiltonian in terms of a Fermi-Hubbard model subjected 
to an artificial SU(2) field llSTl . corresponding to tunneling matrices with position-dependent 
overall hopping rate. 

We will not consider here the fermion doubling problem as it is inessential for the 
discussion of gravitational effects |H Feasible solutions of such a problem in OL simulations 
of 3+1 fermions have been very recently proposed in [66] and BP . 

The motivation for simulating the Dirac equation in curved space-time is at least two- 
fold. The propagation of fermions in curved space is not experimentally accessible in high- 
energy physics/cosmology as the (piece of the) Universe that we are nowadays able to probe is 
practically flat ll74"l , The lattice realization proposed in the paper is an appropriate description 
of the continuum physics for processes occurring over scales of several lattice spaces. The 
realization with cold atoms in OL provides us with "lap-top experiments", paving the way 
to observe exotic effects like the Thermalization Theorem, also known as Fulling-Davies- 
Unruh effect ll75l l76l 1771 (for a review on the subject see ll78l ). Roughly speaking, it states 
that an accelerated observer perceives the Minkowski vacuum as a thermal bath. This is a 
manifestation as Hawking radiation |f79l of the same phenomenon: namely the existence of a 
non-trivial Bogoliubov transformation between the Minkowski vacuum and a space-time with 
a horizon. The effect of the latter, in Quantum Information language, is to "trace-out" part 

\ The mixture of the two Dirac points due to the gravitational background potentially induces a gauge field coupling 
to the composite fermions, (for graphene like lattice see for example I71II421I72II73I ). However, the field contribution 
is relevant in the presence of conical singularities, disclination or dislocation in the graphene language, while is 
negligible when the metric is smooth. 
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of the system, giving raise to a thermal reduced density matrix. Another exotic gravitational 
effect to observe might be the curved space-time version of Zitterbewegung [801 , to compare 
with its recent observation achieved in flat space (54). 

On other hand, simulating the Dirac equation in curved spacetime gives rize to the 
possibility of modelling the analogues of graphene ripples on a flat, square lattice. Such a 
possibility could help in disentangling the action of ripples (which admit a natural "geometric" 
interpretations [8Tj or a gauge field interpretations, e.g. [82]) on the carrier density from other 
contributions. 

Different approaches to quantum simulation include the exact simulation of the dynamics 
of a strongly correlated systems by unitary gates [83), or the the creation of interesting 
strongly correlated states which are ground states of interesting Hamiltonians 11841 . 

The paper is organized as follows. In section [2] we demonstrate how the lattice Dirac 
Hamiltonian in flat space-time 11831 can be obtained as the discretization of the continuum 
Dirac Hamiltonian H. This Hamiltonian is identified then as the Hamiltonian of the SU(2) 
Fermi-Hubbard model considered for instance in [|6T|, where the hopping terms are given by 
the Pauli matrices Cj, i = x,y times a constant hopping rate J. Although the discretization can 
be done at this stage in many different ways, we obtain it by coarse graining the "symmetric" 
formulation of H in the fermion field if) and its conjugate ip\ This strategy turns out to 
be very convenient in order to find the lattice version of the Dirac Hamiltonian in curved 
space-time. As the first exercise, after reviewing the form of the continuum Hamiltonian on a 
generic manifold (admitting a "time" isometry) in section[3] we successfully apply the above 
strategy to the interesting case of the 2+1 Rindler Universe 11861 in section [4] The resulting 
lattice Hamiltonian differs from the Hubbard Hamiltonian of the flat case in the following: 
the hopping rates exhibit a linear dependence with position. Such a simple form is due to 
the cancellation of the spin-connection, once ip and ijj' are treated in the same manner: as a 
consequence the hopping matrices need not be locally rotated. We find that this simplification 
happens in fact for any static spacetime. For the details of the derivation we refer the reader to 
the |Appendix B| We discuss under which conditions, and how the fermion propagation in such 
space-times can be engineered and detected on optical lattices in sections [5] and [6] that are the 
heart of the work. There, several different experimental ways of implementing the Hubbard 
model of interest are presented. We propose the density of states as the relevant observable to 
capture the Dirac physics. It is a measurable quantity both in graphene [87, 88, 89 90 1 and 
in OL (9H [92] |93] |9l] |95] |96] [97] EH experiments. We compute the theoretical value of the 
density of states analytically in the continuum limit using perturbation theory. 

We collect our concluding remarks and discuss further developments in section [7] 

2. Dirac's Hamiltonian on a Lattice 

In this section we show, by discretization the spacial coordinates, that Dirac's Hamiltonian in 
2 + 1 dimensions is a hopping Hamiltonian with non-Abelian tunnelling matrices. To this end, 
let us recall that Dirac's equation in 2+1 dimensions reads j a d a ip = 0, a = 0, 1, 2, where j a 
are Dirac's matrices satisfying {j a , 7&} = 2ry a b, r\ a \> is the metric, and ip is a two component 
spinor. The time evolution of the field ip is, 




(1) 
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from which we can easily read-off the Hamiltonian. As was first pointed out in (85), in 2+1 
dimensions and on a lattice of spacing A the discretized version of equation ([1) is 

(7 )" 1 

id l/j = -i 2 ^ (7 1 {lp m +l,n ~ i>m-l,n) + (V'm.n+l ~ VVn.n-l)) , (2) 

with i/j(x) = ip(mA,nA) = ip m . n - Rewriting the Hamiltonians H = J dxdy ip''H'ip we 
obtain on the lattice, 

H = -i 2A (^i+MlVm,™ + ^m, n+ l7 : Vm,7i) + h - c - ( 3 ) 

m.n 

It is easy to see that ® is nothing more than a Fermi-Hubbard Hamiltonian with non-Abelian 
hopping matrices where the interactions and the effect of the trap have been neglected. In the 
notation of ll6T1l we have U x — 17° 1 j 1 = ia\ and U y — 27° 1 7 2 = ioi (which implies 
70 = 203) and the Abelian flux <E> = 0. 

For further reference, we note that the lattice Hamiltonian of (0 can be alternatively 
obtained by computing 

H = J dxdy ip*Hip = ^ J dxdy (Hip) ] 'if) + ^ J dxdy ^Utp , (4) 

with the substitution of spacial derivatives with finite differences over one lattice space A, 

dA _> fa^M and dy1 p Vw+^-Vw , 



3. Dirac's Equation in Curved Spacetime 

In this section we review the formulation of Dirac's equation on a curved spacetime. Let M 
be an arbitrary curved manifold and let us define a vielbein e", or set of vectors that form 
a basis of the tangent space Tm at each point of M. Here the index \i labels the spacetime 
component and a simply labels the basis vector. 

Although e° is not constant in general, we require it to be covariantly constant (see for 
example l99l ). We introduce a connection u> such that 

This defines u> and the covariant derivative D^. Dirac's equation for a Fermi field ip will be 

ro^ = , (6) 

where 7^ are the curved spacetime gamma matrices. The flat-space gamma matrices j a 
and the 7^ are related by 7^(35) = e£(x) 7a . Thus,we have { 7/1 (a;), 7l/ (a;)} = 2g llv (x) if 
{7a, 7b} = ^f]ab and the vielbein forms an orthonormal basis e°e hM = rf b , with g M „ and rj ab 
the curved and flat-space metrics, respectively. 

By separating the time component we have, j^D^ = ^dt+j^uj^jab+Ydi+j^^lab- 
Now, we are ready to write out the time variation of i/j, 

id^ = -i ( 7 4 ) - 1 Udi + \^f lab + \^t lab \ il> , (7) 

H 

where -f ab = -f a j b . 

Hence, we have identified the Hamiltonian for a Dirac fermion on a curved manifold. In 
section [5] we will proceed to its diagonalization analogously to the flat case and discuss its 
implementation on an optical lattice. 
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4. An example: Dirac Hamiltonian in Rindler's Universe 

As an example, we will derive Dirac's Hamiltonian in a particularly simple spacetime. Let us 
start by recalling that the Hamiltonian in curved spacetime, when a timelike killing vector 
is present, is the integral of the Hamiltonian density on a timelike hypersurface. In 2+1 
dimensions we have, 

H = J dT, 2 ^*Hip . (8) 

The differential is the volume element on a timelike slice and includes the determinant of the 
metric, dT, 2 = ^f^gdxdy. 

In Rindler's space, the metric takes the form 

ds 2 = -{axfdt 2 + dx 2 + dy 2 . (9) 

We are interested in this metric for two reasons. First, it is the metric seen by an observer 
in constant acceleration in fiat spacetime. Therefore, it could have implications for earth- 
dwelling detectors observing cosmic background radiation. Second, it is the near-horizon 
metric of a Schwarzschild black hole. 

The Rindler metric suggests the choice of dreivein e° = |ax|dt, e 1 = dx, e 2 = dy. 
Using © we may compute the spin connection, whose only non-vanishing component is 
in" 1 = a-r-r- Dirac equation (jTJ, greatly simplifies in this spacetime, 



ir)tih = — in\x\ I — 'Vo I rL 4- 

2\x 



id t tp = -ia\x\ ( -72 ( d x + —r ) + 71 9j, ) ip ■ (10) 



dxdyx({d x ^)a x i(j + (d y ip^)a y t/j - ^a x d x ijj - ip^a y d y ip) . (11) 



H 

where we have used that 707172 = — 1, which holds only in 2+1 dimensions. In what follows, 
we adopt the gamma matrices representation choice 

<J X = -72, cr. y =71, a z = -170. 

In order to carry out the discretization analogously to how it is done for the case of a gauge 
theory, we note that the Hamiltonian can be written in terms of the Hamiltonian density dTOb 

as 

H = j dxdy ^ Hip = - j dxdy{H^)H + \ j dxdy ^ Hip 
ia 
~ ~2 

In this symmetric form the spin-connection term disappears and it turns out again once 
integrating by parts. 

The discretized version of the Hamiltonian (fTTl i is simply obtained by the substitution 

d x 1p — > + ^' m - n an( J Q yi p — » V'm.n+^-^m.,! ^ w J tn x _ ^ m an( J y _ rea( JjJy 

gets 

H= J^2 m {^In+l^^ra.n + 1p m ,n+lVy1p m ,n^ + H.C. . (12) 
m,n 

Therefore, a lattice with hopping matrices given by JjZT \ = iraui, i = x,y, growing linearly 
in the x direction gives an appropriate description of free massless fermion in a Rindler 
spacetime. Such an Hamiltonian can be in principle implemented in a OL. This problem 
will be tackled in the next section. 
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5. Dirac equation in curved spaces and optical lattices 

In view of the explicit realization in a OL we focus to the classes of spacetime where the 
massless fermion propagation can be described by the Hubbard model of the form 

= 7 V J mn ypL+l^^m^n + ^n+l^Vw) + h.c. . (13) 



2 ^ 



As discussed in full details in the Appendix B we find that the lattice Hamiltonian reduces to 



this simple form only if the metric is static. In 2+1 dimensions it is equivalent to say that it 
can be written (in a certain coordinate system) in a diagonal form as 

ds 2 = -f(x,y)dt 2 + f- 2 (x,y)e 2 ^>y\dx 2 + dy 2 ) . (14) 

With this choice of coordinates the hopping rate is simply given by J mn = e®( Xm ' Vn \ 

This is not the only requirement to be satisfied to reproduce such propagation on a OL, 
however. Indeed, it is important to note that, even if the Dirac Hamiltonian written in the 
symmetric form ( fT3l is the same for any choice of the function / in ( TBI ), the corresponding 
Hamiltonian system is distinct for each metric as the canonical momentum is 

H = iV^gi'J* = if^e 2 ®^, 

and depends explicitly on /. Now, in optical lattice experiments the canonical momentum is 
fixed by the anticommutation relation to be simply iijp that implies / = e . Thus the cold 
fermions in the optical lattice simulate the propagation of massless fermions in a metric of the 
form[§| 

ds 2 = -e^dt 2 + dx 2 + dy 2 . (15) 

The Rindler metric © corresponds to $ = log(aa;). 

Before moving to the explicit implementation of the Hamiltonian (fT3l l in OL, let us 
briefly discuss to what extent it is a good description of the dynamics of massless fermions 
in a space time given by (I15t . There are two kind of limitations. On one hand, due to finite 
size of the OL we are able to cover only a finite portion of spacetime. On the other end, the 
lattice approximation is valid when the metric, or the function <£>, is sufficiently smooth and 
slowly varying over one lattice space A. It is worth to note that the second problem can be 
circumvented by using techniques from lattice gauge theory (the continuum limit is obtained 
by extrapolation). 

These two limitations should not obscure the physical content, however, as OL with up 
to 300 x 300 sites can be achieved. This implies that the overall variation of the metric over 
lattice can be of order one. Further considerations in the case of Rindler space are given in 
section|7] 



5.7. Experimental realization 

In this section we discuss briefly how one can achieve the appropriate site-dependent hopping 
rate J mn in an OL. We propose two different techniques. 

Recently it was proposed in 11 1 II [311 . that the intensity of the hopping can be 
tailored almost at will by considering bichromatic spin-independent superlattices that trap 
the hyperfine states of alkali bosonic or fermionic atoms. The split of the hyperfine levels is 
controlled by a magnetic field. The hopping between neighboring Zeeman sublevels of the 
F (lower) hyperfine manifold, i.e. our "electrons", is induced via adiabatic elimination of an 

§ It is worth noting that for a generic function <E> the metric is curved and not Weyl-invariant, cfr 11001 . 
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intermediate F = F ± 1 (upper) hyperfine manifold coupled to F via an off -resonant Raman 
transition. For the details of the scheme for fermionic 40 K atoms, we refer to the original 
proposal J3TJ . 40 K have F =9/2, and allow in principle to simulate any pseudo-spin F 1 <9/2, 
employing the splitting of the Zeeman sub-levels in a magnetic filed, and optical pumping to 
the relevant sub-levels. For the purposes of the present paper it sufficient to have F' = 1/2, or 
alternatively to use from the very beginning atoms with F = 1/2 in the ground state manifold, 
such as 6 Li. 

To be more concrete, consider a 6 Li Fermi gas loaded on a 2D square optical lattice of 
size LxL, where the relevant information of our quantum simulator is encoded in the Zeeman 
sub-levels of the hyperfine manifold F =1/2. Laser-assisted tunneling methods allow us to 
design arbitrary operators U ru dressing the hopping between lattice sites r — > r + u, where 
r = mx + ny, m,n,E {1...L}, and v 6 {x,y}. Usually, such schemes rely on Raman 
couplings to auxiliary states trapped in the links of the original lattice and belonging to a 
different hyperfine manifold. Here, following IU01II3T1 we use bichromatic spin-independent 
superlattices , and use the secondary minima of F =3/2 as bus states to mediate the hopping. 
Note that the individual addressing of each hopping rate is granted by the Zeeman splitting 
within the hyperfine manifolds, and the different detunings of the Raman lasers. These 
detunings can be quite large, so that the lifetimes of atoms on the lattice (limited by photon 
absorption and spontaneous emission) can be quite large, of order Tj ~ Is. By making the 
Raman laser intensity/detunings and/or Zeeman splitting spatially dependent one obtains the 
desired spatially dependent hopping rates which is necessary for the realization of curved 
space-times, equation d 1 3t . On top of that, is is possible to use Feschbach resonances to turn 
off the atom-atom scattering, and make the system essentially non-interacting. 

An alternative, and perhaps even a simpler method to realize the Fermi-Hubbard model 
of the form of ( fT3l can be achieved by taking into account the finite waist of the lasers 
used for the generation of the hopping terms. In general, this is an undesirable feature, and 
it can usually be neglected. Indeed, typically Gaussian laser beams of waist w are used, 
characterized by a Rayleigh length — the the distance along the direction of propagation from 
the waist to the place where the area of the cross section is doubled — zr = irw 2 /\ with A 
denoting the wavelenght. That is to say, within a volume of w 2 x z r , the ideal planar wave 
is a good approximation, at least around the center of the beam. For a lattice with L = 30, 
its linear size is 30 x A/2, so that focusing a laser on the whole lattice w ~ 15 A leads 
to zr ~ 700A, so that the plane wave description can be even used for an array of a few 
hundreds of 2D 30 x 30 lattices. 

There are, however, no technical obstacles to focus the lasers on much smaller spots, 
smaller than LA/2, keeping z r still quite large. Actually, such spacial modulation of the 
intensity due to the waist of "real" lasers can be used to induce hopping terms that depend 
non-trivially on the position (29l . In general, the hopping rate, i.e. the modulus of the 
hopping term, is proportional to the intensity of the laser producing it. For instance, taking the 
paradigmatic example of the hopping induced by the Raman transition in Jaksch and Zoller's 
setup 1 131, an optical lattice implementing the Hamiltonian of dT3l > can be achieved for Raman 
lasers propagating all in the same direction, once we consider only the radial waist and neglect 
the waist along the beam[jj|. In this case, the shape of the laser intensity will correspond to the 
e* factor of the metric. 

Comparing the two methods explained above, the former has the advantage that, in 

|| At the moment, it is not clear to us whether it is possible to engineer the hopping matrices cr x and cr y with parallel 
and anti-parallel Raman laser. Actually, we do not know of any explicit realization of such hopping in a setup a 
la Jaksch and Zoller. Maybe, the dark state or slow-light method |14| is more promising. The main point that the 
hopping rate is controlled by the intensity, remains valid. 
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principle, any shape of J mn , i.e. any metric of the form ( fT51 ), can be engineered, but at 
the price of dealing with a quite involved experimental apparatus, while the latter, although it 
allows for a restricted choice of J mn (for instance a Gaussian shape) is almost for free. The 
desired hopping rate profile is obtained by reverse engineering of laser waist. 



6. Density of states at the Fermi level as simple observable 

Let us turn to the discussion of possible detection schemes of Dirac physics in curved 
spacetime. A simple observable that characterizes such physics, and that contain information 
about the effects of the beam waist is the density of states RTl 11021 11031 . This quantity 
is routinely measured in graphene using scanning tunnel [89 90 1 and electron transmission 
spectroscopy [87, 88 1. For ideal graphene at zero temperature the density of states is zero at 
the Fermi level and is proportional to \E\ (once we have fixed E F = 0) as the charge carries 
are described by massless fermions propagating in fiat space. In the presence of deformations 
of the graphene sheet its deviation from the free behavior can be analytically computed at 
first order from the propagator of the Dirac equation. Following 1811 . by modeling such 
deformations as perturbations of the Minkowski metric it is possible to compute the Green 
function treating the correction to the free equation as an interacting term V. In fact, in this 
section we will first reproduce the computation of lIHTl for a metric of the form (fT5l l, instead 
of the spatial deformation (as in (IA.1I )) considered there. 

Our final goal is to determine the local density of states, defined by 



dw^/— gd rp(w, r) = # of states , (16) 
in terms of the Feynman propagator (see appendix) using the relation 

(17) 



p(w, r) = sign(iu)— Im Tr S F (w, r, r)7* 

7T L 



where r = (x, y). 

In order to find S F (w, r,r) we start by using the defining equation for the fermion 
propagator 

-iV 3 57 M ^r^F(x,x')-^ 3 (x-x'), (18) 

where the x indicates a point of the spacetime, to a metric of the form ( flST l. By retaining terms 
linear in <1>, the above equation can be written as 

- (i(7 M ^r)f fat + V(x))S F (x, x') = 5 3 (x - x') , (19) 



where 



V = l 1 l[ <S>d X + -d X Q ) + i 72 ( §dy + -By® ) , (20) 



is the effective "external potential". As we are interested to S F (w, r, r') and due to the time 
translation invariance 

S F (x,x') = S F (t-t', r,r') = J |V™('->')£ F Kry), 

it is convenient to perform the Fourier transformation in time of $1% 

(w l0 - *V r ■ 7 - V(r)) S F (w, r, r') = <5 2 (r - r') . (21) 
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The above equation can be solved consistently within the first order approximation by 

^Kr,r') = J A 2 r"S%{w,vy')V{v")S%{wy,v'), (22) 



where 



J iZttY ur — + te 



(23) 



is the free fermion propagator and space translation's invariance holds. By using the Fourier 
transformation of <£>(r") 

d 2 p 



$(r") 



(27T)' 



;e lp r *(p) ; 



equation (l22l can be explicitly computed performing the integration in r". 

The relevant contribution to the trace turns out to be linear in w. Explicitly 

}2 T 



Tr[^Kr,r) 7 „] = J AjLj*>-'*(p)r(w, p) , 



with 



i>, P ) 



d 2 & 



4tu|k- ip| 2 



(24) 



(25) 



(26) 



(2tt) 2 (w 2 - k 2 + ie)(u; 2 - (k - p) 2 + ic) 

The above integral is logarithmically divergent, but its imaginary part is not. It is easy to 
show that this is the only part contributing to the density. Indeed, as <&(p) is the Fourier 
transformation of a real function, < E > (p)* = <&(— p), and r(to,p) is even in p, T(w, p) = 
T(w, — p), one finds 

d2 P p -«p- 



(Tr[^(u;,r,r)7 



which immediately implies 



5p(w) = sign(iu)- 



d 2 r 



r <I>(p)*r( W ,p)* 
( 27r ) 

° e ipr $(p)Imr(w,p) 



The explicit expression for Im T(w, p) is 



w 



Imr(w, p) — 2w Im 



(2^) 2 
arctanx(w,p) 



(27) 



(28) 



(29) 



X(w,p) 

, and p = |p|. More details are given 



where y(w,p) = -, p 

V ' , J\4w 2 -p 2 \(e(iw 2 -p 2 )-iQ(-<iw 2 +p 2 )) 

in the | Appendix D| 

Hence, the density of states always receives a correction proportional to $(r) itself 
r . 9.U,>I Unl r . r r> 2 

Sp(w 



y d vM* (r )-My d v/A^e-*(p)Im 



arctanx(if,p) 



(30) 
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6.1. Example: Density of states for a Gaussian beam 

We are now able to compute the correction to the density of state when the hopping J has 
a Gaussian shape due to the finite laser beam waist. Under the assumption that the Raman 
lasers propagate along y direction, this implies that $>(x, y) = —(f) 2 , where a is of order 10 2 

lattice spacings fjp. As <&(p) = + ^2 S" {p x )8{p y ), the second term of d30t turns out to be 
zero and the correction to the density of state simply reduces to 

6p(w) = J d 2 r^U(r) (31) 

providing a clear experimental signature. The local density of states gets a quadratic 
correction in x. 

Incidentally, the same cancellation happens in case of an exponential behavior of J, i.e. 
for a linear 3>(r). Consequently, the simple relation (T3TT > applies also to this case. 



6.2. Experimental detection 

A recent review of the detection methods that can be applied to investigate Dirac physics 
with ultracold fermions in non-Abelian gauge fields is contained an article authored by one 
of us |9T| . Here we just summarize this discussion with particular focus on density of states. 
Let us start by observing that for a non-interacting Fermi gas at T = 0, the total number of 
fermions Np = J' 1 dEp(E), where p, is the chemical potential, equal at T = to the Fermi 
energy Ep. We see that piEp) = dNp/dEp so that measuring of the variance of Np with 
Ep allows to determine p(Ep). If the systems is confined additionally in a slowly varying 
harmonic potential V(r), a local chemical potential can be introduced /i(r) = Ep — V(r), 
and the corresponding local density of states, related to the local density by 

n(r) = J dEn(E)Q(p(r) - E), 

where 0(.) is the Heaviside (step) function. In this case we get 

dn(r)/dp(r) = p(p(r)), 

i.e. a similar formula to the Streda formula used in ||9l1 for detection of Hall conductivity. 
The determination of density of states can be achieved by: 

• Measurements of the total number of fermions as a function of the chemical potential. 
Here, the best currently available methods are: direct in situ individual atom detection 
||92][93][94], or quantum spin polarization spectroscopy |95ll96l 

• Measurements of the (coarse-grained) local density of fermions as a function of the 
local chemical potential. Again, the best currently available methods are: direct in situ 
individual atom detection |92ll93ll94l , or quantum spin polarization spectroscopy with 
spatial resolution |97l . 

• Measurements of frequency-momentum resolved single particle excitation spectrum, 
such as those being done in Bragg (Raman) scattering spectroscopy (for a state-of-the-art 
report see j98l ). The spectrum in such processes is proportional to the density of initial 
states of the scattering process. 

Of course, many other methods, such as atomic ARPES, noise interferometry, or even 
absorbtion and/or phase contrast imaging can give at least indirect information about p(E). 
All of these methods are well developed in experiment with ultracold atoms (see ||9T1 and 
references therein). 
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7. Conclusions and Outlook 

In this paper, we discussed the simulation of the Dirac equation in artificial curved space- 
time with cold atoms. We showed that using state-of-the-art techniques it is possible to 
simulate relativistic fermion dynamics in curved spacetimes with a. flat 2D square lattice for 
an interesting class of 2+1 metrics. Moreover, we pointed out the relation between a certain 
class of Hubbard models and Dirac's Hamiltonian in curved backgrounds, which can be 
employed to make analytic computations in the continuum limit of the former. We proposed 
to characterize the Nature of Dirac fermions on the lattice by measuring the density of states 
at the Fermi level. This observable can be, on one hand, analytically computed in perturbation 
theory in terms of Dirac propagator, and, on the other hand, is accessible to measurements. 

The present study opens the way to the direct observation of elusive effects such as 
Rindler noise. Because we deal with odd dimensional (2+1) Rindler system, the Dirac 
thermal noise, measured by an ideal point-like De-Witt detector as a consequence of the 
local acceleration, is expected to be "anomalous" (see Ch.8 of l78l ). i.e. it should follow 
Bose-Einstein distribution. This issue is currently under investigation. 
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Appendix A. Dirac Hamiltonian in spatially (graphene like) deformed metric 

Now we consider a different situation where the 2+1 metric is spatially non-trivial. Such a 
case is relevant in describing the properties of a graphene sheet with ripples. The most generic 
spacial deformation (at least in some patch) can be always written as 

ds 2 = ~dt 2 + e 2 ^y\dx 2 + dy 2 ). (A.l) 

The driebein are 

e° = dt e 1 = e^'^dx e 2 = e^'^dy i=x,y, 

and the spin-connection can be chosen to be non-trivial in the spacial part only 

u; 12 = dy^dx - d x 1>dy. (A.2) 
It follows that the curvature is 

Q = 2e- 2 *n 12 y = -2e~ 2 *(d 2 + 9 2 )$ . (A.3) 

For instance, the slices of the metric ( IA. 1 b at constant time will be spheres or hyperboloids 
for $ a positive or a negative quadratic form of x and y, respectively. 
Applying (0 to this case we find 

iditjj = Hi' = -i^e"* (ji{d x + \d y <S> ll2 ) + y 2 (d y - ^d x ^ 12 ) j j> . (A.4) 
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We are tempted to interpret the above Hamiltonian as that of fermions coupled to a 
"geometric" non-Abelian vector potential A = (d y &, —d x &)o~ z . Adopting the the gamma 
matrices' representation of the previous section, we can write it as 

nijj = -ie~* (a x (d x + iA x ) + a y {d v + iA v )) if> , (A.5) 

where the presence of A indicates that the rotation in the xy-plane, i.e. the 5*0(2) subgroup 
of the Lorentz group, is promoted to a local symmetry in the background described by the 
metric dA.U . Such identification is related to the treatment of the conical defects in the 
graphene sheets (dislocations and disclinations) as sources of magnetic fluxes. However, that 
this interpretation here is misleading as the gauge group does not commute with the spacetime 
symmetry, as a z anticommutes with a x and cr y . 

Taking in account this fact, the Hamiltonian density can more appropriately be written 

as 

H (s . d .) = -i (a x (d x + ^d x $) + a y (d y + ^d y $)j V , (A.6) 

where the symmetric role of x and y is evident. 

Now we are ready to compute the total Hamiltonian. By rewritting (0 we get 

H= [ Axdye 2 *^^U (s . d ^. (A.7) 



Using the the same manipulations as in the previous section we can recast it into a form where 
the spin-connection is not present, 

H =\\ dxdye^'V' 1 ^ ({d,^)^ - ^ a^) . (A.8) 

i— x,y 

It follows that the discretized version of H, as in Rindler spacetime, takes the form of a SU(2) 
Fermi-Hubbard model with the modulus of hopping depending on the position 

i v-^ e #(mA,nA) , . 
IFH = 2 A yPm+l^x^m^ + 4>m,n+l<?ylp m ,n) + h.C . (A.9) 



Hp 



At first sight it seems very surprising that the discretized Hamiltonian of massless fermions 
in Rindler geometry coincides with the one of fermions propagating in metric of the form 
of ( I A. It . Indeed, by taking $ = ln(aa;) the expression JA.9t reduces to ( fT2b . Such an 
apparent contradiction disappears upon closer inspection. At the end, from the point of view 
of the Dirac Hamiltonian for both metrics, what has changed with respect to the flat case is 
the effective speed of light, or equivalently the hopping rate, which becomes position (and 
direction) dependent. 

Although the origin of the position-dependent hopping rate is different in the two cases, it 
comes from the Hamiltonian density in the Rindler case while is due to the invariant measure 
in the other, the effect is the same. Roughly speaking, there are two possible way of modifying 
the effective speed of light in one spacetime direction, let us say x: one is to change g u while 
the other is to change g xx by an inverse factor. 

Nevertheless, the eigenfunctions and the spectra of the two Schrtidinger problems remain 
different as the Hamiltonian densities in the two cases are. 
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In order to treat the generic case let us retrace a few steps and analyze the formal expression for 
the Hamiltonian (|8}. First of all, we show that it corresponds to the Legendre transformation 
of the relativistic Lagrangian: 

l= fc = i[ V^di'^D^, (B.i) 



once we have chosen the coordinates to have that the timelike Killing vector is K = dt, which 
implies d^g^ = 0. Indeed, by defining % = j^pdtTp — the expression ([8]l is recovered 

1 



H = Jh = -*JV-9i> [l l d l + j7 M V 6 7a6 ) V> , * = x, y. (B.2) 

Now, it is instructive to check explicitly that the Hamiltonian above is a Hermitian 
operator due to the existence of a timelike isometry. By using that (7'*)' = 7o7 M 7o and 
noticing that 

$ 7" W ^ = W W , 

we get that 

1 



Hi = i j V=9 - w^labl^J V> ■ (B.3) 

In order to compare the above expression with H it is convenient to integrate by parts and 
rewrite it as 

# f = - ij y/=g$ (l% + \w ll ab 'Yab'V M + 1% In J~g + drfj Tp 

= H-i J V^i'Qv b [7a^7 AI ]+7 l 9 i lnV^+9a 1 )^- (B.4) 

W and H are the same if and only if the metric, and consequently the driebein, are 
invariant under time translation. Indeed, due to the conventional constraint V M e° + iy„ a e£ = 
0, it is true that 

e^eg = -egV^e" = -e va w^ a e" b , (B.5) 
which, using the commutator [y a b, 7^] = 4e AIC 7[ a ?7b] c , implies 



(B.6) 



The cancellation follows from the identity 

d a l + 7*$ ^ sl~~g = V l7 l = V m t" , (B.7) 

as <9 M In = T" and d t 7 * = j*dt In \J~^g = if and only if d t g^ v = 0. 

At this point, we can use expression ( IB. 3b for W to get a Hamiltonian symmetrical in t/j 
and ip^ . By writing H = + iff) we find 

H= l -j^—g UdiWI) - fa'drf) - ^V 6 {7oi, 7"}V>J ■ (B.8) 

Let us characterize the term 8 = — ji^w^ilab, 7^}, which can be regarded as the 
obstruction to write the lattice Hamiltonian simply as 

H = \ ^2 Jmn (y > ti+l,n <T z'<Pm,n + V'm.n+l^l/V'm.n) + h - c - > ( B - 9 ) 



Dirac Equation For Cold Atoms In Artificial Curved Spacetimes 



15 



as found for the special cases discussed in the previous sections. Using the relation between 
the spin-connection and the driebein and the anticommutator {j a b, 7 M } = 2e^ c 7 a 6 c -in 2+1 
dimensions it reduces to {-f a b, 7 M } = — fe^e^c, it follows that 

6 = ^ a e vh d^ v e abc = ^e ia e vh diel e abc . (B.10) 

Hence, we conclude that 6 is identically zero when the metric is diagonal, as for the Rindler 
metric (0 and the spatially deformed metric ( IA.U . In order to be as general as possible, we 
observe that any 2+1 spacetime admitting a timelike Killing vector can be always reduced, at 
least locally, to the form 

ds 2 = -f(x,y)(dt + X(x,y)) 2 + f- 2 (x, y)e 2 ^\dx 2 + dy 2 ) , (B.ll) 

where X(x,y) = X x (x,y)dx + X y (x,y)dy is one-form independent of time. 
It follows that: 

e ° = f{dt + X) , e' = f~ 1 e^dxi , xi = x,y, i = 1, 2 , 
and that the inverse driebein are 

eo = f~ X dt , ei = -fe~®X x d t + fe~*d x , e 2 = -fe~^X y d t + fe~*d y 
The spin-connection can be chosen to be: 
w oi = fe -* dxfdt + tl_ (d x (fX y ) - d y (fX x )) dy 

W 02 = fe-^dyfdt - (d X (fXy) - dy{fX X j) dx 



2 



+ (dy($ - In/) - X x (d x (fX y ) - d y (fX x ))^j dx 



d x ($ - In/) + X y (d x (fX y ) - d y (fX x )fj dy. (B.12) 

This implies the following relation for the coefficients: 

01 02 J--2 * 12 

w y = -w x = f e w t , 

and that w y 12 can be obtained from w x 12 by exchanging x with y and 1 with 2 (which amounts 
for the minus sign), accordingly to the symmetry of the metric ( IB. 1 Q . 
We are ready to compute O. After some algebra we get 

f 3 -2* 

9 = J —^— (d X Xy - dyX X ) . (B.13) 

The above expression means that the spin-connection will cancel in the symmetric form of 
the Hamiltonian if and only the metric is static, i.e. the off-diagonal terms due to A can be 
reabsorbed by a change of coordinates. To prove this we note that the condition 8 = implies 
that the form A is exact, i.e. there exists a function F = F(x, y) such that dF = X. Indeed, 
this is the case as as can be seen by redefining the time coordinate as T = t + F, as the final 
metric is diagonal. 

With the choice of the metric (IB. lib , the overall hopping rate J mn is 

,J mn — c 
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Appendix C. The density of states 

This section is devoted to deriving the relation between the density of states as a function of 
the energy p(E) and the propagator in a spacetime of dimensions (1, d). In order to do so we 
will follow a constructive procedure. By definition we have 

J &E p{E) = # of eigenstates of the Hamiltonian H . (C.l) 

We will express the number of states as the number of poles in the propagator when averaged 
over the set of eigenstates of H, {|n}}, which we demand be complete and normalizable (as 
is the case for any sound Hamiltonian operator). Furthermore, we assume that the position 
operator eigenstates |r) can be completed to give an orthonormal basis |r)|£) where i encodes 
all internal degrees of freedom, such as the spin. Using the theorem of residues and integrating 
on a rectangular contour around the real axis of width 2e > , we can write: 

-^UW l l — 

-2m J ^\E + ie-E n E-ie-E 
1 , , 



#e.s. of H = 



ilm / d E J2( 



E — H 



ilm [dE [ d d rdV g (n|lr)|«)(r|(«| g _ ^ . e \r')\i')(v'\(i'\\n) 



7T 



!)././' 



= - -lm f dE [ d d r{v\(j,\ - )— lr)K) . (C.2) 
7T J J E - H + le 

Form the above equation we conclude that 

P {E) = ilm / d d r{v\{i\ 1 . \r)\i) , (C.3) 

7T J H — E — It 

that is the equation 34 given in 11 1 041 . To warm up we compute the density of state in the free 
case. As the Hamiltonian is diagonal in momentum space, it is convenient to write p{E) in 
this basis: 

P (E) = Iim Jd 2 n p (<p,+I^t-Ip>+> + ^e-JTTT^- 

= > J* 2 rJ^ p ((P.+lk) |p|2 _]g K f ign(g)£ (r||p ; +) 

. .. . -Ipl +E . .. 

+ <P> ~H r /]pl2 ^ -T7T77^-( r IIP' 

1 T f 2 f d 2 p 2E 



7T 



p| 2 - E 2 - isign(E)e 
Im / d 2 r / 77^ , (C4) 



(2tt) 2 |p| 2 - E 2 - isign(E)e ' 



where d2^l p = j2wT*2\p\ ^ s Lorentz invariant measure and the normalization of the 
momentum states |p, ±) is fixed accordingly. By taking to be S — J d 2 r the volume of 
the system and by going to polar coordinates one finds 

. , S f +co dp 2pE 
p{E) = -Im / -f- y ,. . (C.5) 

7r Jo 2tt p 2 - E 2 — isign(E)e 

The last integral can be solved in many ways, for example by changing variable to z = p 2 and 
regularizing the integral with cut-off 
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f A dz E _ E f A 2 -E 2 ( 1 - zsign(£)e \\ 

J 2-kz-E 2 - isign(E)e ~ 2tt ^ ° g( E 2 ' + ° S \ 1 + isign(E)e ) ) 

= fr ( log (^-^") + 2 ™( (£) - ®(-^))) CC-6) 
As the imaginary part is independent of the cut-off, the final result, as expected, is 

p(E) = S 1 -^ , (C.7) 

7T 

which is in agreement with the result quoted in 11041 . 

It is worth-while to note that p(E) can be written in terms of the Feynman propagator. 
This fact can be derived in a more general setting. By definition, E _^j +ie is the Fourier 
transformation in time of the retarded propagator defined by the equation 

(id t -H)Q + {t-t') = 6(t - t') g> I spin (C.8) 

with boundary conditions Q+(t — t') = for t < t' . On the right-hand side of ( IC.8I ), the 
identity in spin-space has been written explicitly to remind the reader that Q + (t — t' ) in general 
acts as a matrix on the internal degrees of freedom. After multiplying on the left by 7* and on 
the right by 7 t and taking the expectation value with eigenstates of the position operator the 
above equation gives 

- i^d^G+it - t', r, r')7t = S(t - t')6 d (r - r') <g> I spin , (C.9) 

where we use 7*74 = 1. Hence, we conclude that Q + (t — t 1 , r, r')7t is related to the Feynman 
propagator as it solves the same equation. The precise relation can be derived by taking the 
boundary conditions into account. This can be explicitly checked by Fourier-transforming to 
momentum space. Indeed, 

1 E - -ftl l Pi Ej t - 7>i 

-It = 7^ 1 10 , ■ ■ — 77^7* ~ 



E-H + ie 1 — |p| 2 + zsign(£7)e ' E 2 - |p| 2 + isign(£T) £ 

1 



,pr- (CIO) 

The above relation between the retarded propagator and Feynman propagator can be extended 
using perturbation theory to the interactive second quantized formalism. To conclude, let us 
remark that it is better to use the local definition of density ( fTTI i because it is easy to make it 
generally covariant in order to apply it in a curved gravitational background. In this way, the 
generalize notion of inner product is properly taken into account due to ( fT9b . 

Appendix D. The computation of T 

In order to compute Imr(w, p) we note that the integral of ( f26T > can be split in 

f d 2 k ( 1 w 2 -|p 2 +p-k \ 

r( w , P ) = 4 W / 7 — - f -a 7 - + — - a - . (D.i) 



(27r) 2 \k 2 — w 2 — ie (k 2 — w 2 — «e)((k — p) 2 — w 2 — ie) 

The first integral has been computed above for the free case and its imaginary part gives 

f d 2 k Aw „ _ ^ s 

Im / - 2 = 2w. (D.2) 
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The second integral is convergent and can be computed using Feynman parameter. After some 
algebra one gets 



Aw Im 



d 2 k 



3„2 



p k 



(2tt) 2 (k 2 - w 2 - ie)((k - p) 2 - w 2 - ie) 



= Im \l Aw 2 + ie — p 2 arctan 

np \ 



\J Aw 2 + ie — p 2 



J\Aw 2 -p 2 \ Q(Aw 2 - p 2 ) Im arctan 

np y 

-Q(-Aw 2 +p 2 ) Re arctan 



\J Aw 2 + ie — p 2 



\J Aw 2 + ie — p 2 



(D.3) 



where \J Aw 2 + ie — p 2 = ^/|4w 2 — p 2 \ (@(Aw 2 — p 2 ) — iQ(—Aw 2 + p 2 )) and p is the 
modulus of p, p = |p|. 



References 

[1] Fetter A 2009 Rev. Mod. Phys. 81, 647. 

[2] Dalibard J et al 2008 Rev. Mod. Phys. 80, 885. 

[3] Lewenstein M et al Advances in Physics 56, 243. 

[4] Goldman N 2009 Quantum transport in lattices subjected to external gauge fields: The quantum Hall effect in 
optical lattices and quantum graphs (VDM Verlag Dr. Muller Aktiengesellschaft & Co. KG). 

[5] Dalibard J et al 2010 Artificial gauge potentials for neutral atoms Preprint arXiv: 1008.5378 

[6] Pitaevskii L and Stringari S 2003 Bose-Einstein condensation (Clarendon Press, Oxford). 

[7] Madison K et al 2000 Phys. Rev. Lett 84, 806. 

[8] Abo-Shaeer J et al 2001 Science 292, 476. 

[9] Cooper N et al 2001 Phys. Rev Lett. 87, 120405. 
[10] Wilkin N and Gunn J 2000 Phys. Rev. Lett. 84, 6. 
[11] Paredes B and Cirac J 2003 Phys. Rev. Lett. 90, 150402. 

[12] Sarajlic E et al 2009 Pushing the frontiers of atomic physics: Proceedings of the xxi international conference 
on atomic physics, edited by R. Cote, PL. Gould, M. Razman, and W.W. Smith, (World Scientific, New 
Jersey). 

[13] Gemelke S C N et al 2010 Rotating Few-body Atomic Systems in the Fractional Quantum Hall Regime 

P»wto larXiv:1007.2677l 
[14] Juzeliunas G and Ohberg P 2004 Phys. Rev. Lett. 93, 033602. 
[15] Jaksch D and Zoller P 2003 New J. Phys. 5, 56. 
[16] Mueller E 2004 Phys. Rev. A 70, 041603. 
[17] Sorensen A et al 2005 Phys. Rev Lett. 94, 086803. 
[18] Tung S 2006 Phys. Rev Lett. 97, 240402. 
[19] Bhat R et al 2006 Phys. Rev. Lett. 96, 060405. 
[20] Klein A and Jaksch D 2009 Europhys. Lett. 85, 13001. 
[21] Ruseckas J et al 2005 Phys. Rev. Lett. 95, 010404. 
[22] Osterloh K et al 2005 Phys. Rev Lett. 95, 010403. 
[23] Lin Y J et al 2009 Phys. Rev. Lett. 102, 130401 . 

[24] Boada O et al 2009 Gauge-Away Effect in Cold Gases on Optical Lattices Prgprinf[arXiv:0909. 39371 
[25] Boada O et al 2010 Simulation of gauge transformations on systems of ultracold atoms Preprint 
larXiv:1009.1552l 

[26] Moller G and Cooper N R 2010 Condensed Groundstates of Frustrated Bose-Hubbard Models Preprint 

larXiv: 1009.44201 
[27] Lin Y J et al 2009 Nature 462, 628. 
[28] Spielman I B 2009 Phys. Rev. A 79, 063613. 
[29] Giinter K J et al 2009 Phys. Rev. A 79, 01 1604. 
[30] Gerbier F and Dalibard J 2010 New J. Phys. 12, 033007. 

[31] Bermudez A et al 2010 Wilson Fermions and Axion Electrodynamics in Optical Lattices Preprint 
larXiv:1004.510ll 

[32] Goldman N et al. 2010 Engineering Time-Reversal Invariant Topological Insulators With Ultra-Cold Atoms 
Preprinf larXiv: 1002.02191 



Dirac Equation For Cold Atoms In Artificial Curved Spacetimes 



19 



[33] Spielman I B Talk at DARPA meeting 2010. 

[34] Wallace P R 1947 Phys. Rev. 71, 622. 

[35] Novoselov K S et al 2005 Nature 438, 197. 

[36] McCann E and Fal'ko V I 2006 Phys. Rev. Lett. 96, 086805. 

[37] Katsnelson tAIetal 2006 Nat. Phys. 2, 620. 

[38] Zhou F and Semenoff G 2006 Phys. Rev. Lett. 97, 1 804 1 1 . 

[39] Geim A K and Novoselov K S 2007 Nat. Mater. 6, 183. 

[40] Novoselov K S et al 2007 Science 315, 1379. 

[41] Matulis A and Peeters F M 2007 Phys. Rev. B 75, 125429. 

[42] Jackiw R and Pi S Y 2007 Phys. Rev. Lett. 98, 266402. 

[43] Pendry J B 2007 Science 315, 1226. 

[44] Cheianov V V et al 2007 Science 315, 1252. 

[45] Castro Neto AWetal 2009 Rev. Mod. Phys. 81, 109. 

[46] Zhu S L et al 2007 Phys. Rev. Lett. 98, 260402. 

[47] Soltan-Panahi P et al 2010 Multi-Component Quantum Gases in Spin-Dependent Hexagonal Lattices Preprint 

larXiv: 1005. 1276] 

[48] Unanyan R G et al 2010 Spinor Slow-Light and Dirac particles with variable mass Preprint arXiv:1005.3429 

[49] Volovik G E 1992 Exotic Properties of Superfluid 3He (World Scientific Pub Co Inc). 

[50] Volovik G E 2003 The Universe in a Helium Droplet (Oxford University Press, USA). 

[51] Bevan TDCc/ al 1997 Nature 386, 689. 

[52] Lamata L et al 2007 Phys. Rev. Lett. 98, 253005. 

[53] Bermudez A 2007 Phys. Rev. A 76, 04180KR). 

[54] Gerritsma R et al 2009 Nature 463, 68. 

[55] Casanova J et al 2010 Phys. Rev. A 82, 020101(R). 

[56] Gerritsma R et al 2010 Quantum simulation of the Klein paradox Preprint arXiv: 1007.3683 

[57] Cirac J I et al 2010 Cold atom simulation of interacting relativistic quantum field theories Preprint 

larXTv: 1006.2975 1 

[58] Juzeliunas G et al 2008 Phys. Rev. A 77, 01 1802. 

[59] Burrello M and Trombettoni A Non-abelian anyons from degenerate Landau levels of ultracold atoms in 

artificial gauge potentials Preprint arXiv: 1006.0453 

[60] Vaishnav J Y and Clark C W 2008 Phys. Rev. Lett. 100, 153002. 

[61] Goldman N et al 2009 Phys. Rev. Lett. 103, 035301. 

[62] Juzeliunas G 2010 Phys. Rev. A 81, 053403. 

[63] Merkl M et al 2010 Phys. Rev. Lett. 104, 073603. 

[64] Ginsparg P H and Wilson K G 1982 Phys. Rev. D 25, 2649. 

[65] Wilczek F 1978 Phys. Rev. Lett. 40, 279. 

[66] Lepori L 2010 (3+1) Massive Dirac Fermions with Ultracold Atoms in Optical Lattices Preprint 

larXiv: 1004.47441 

[67] Lim LKetal 2008 Phys. Rev. Lett. 100, 130402. 

[68] Lim LKetal 2010 Phys. Rev. A 81, 023404. 

[69] Satija Metal 2006 Phys. Rev. Lett. 97, 216401 . 

[70] Montvay I and Miinster G 1997 Quantum Fields on a Lattice (Cambridge University Press, Cambridge). 

[71] Chamon C 2000 Phys. Rev. B 62, 2806. 

[72] Pachos J K et al 2007 Eur. Phys. J. 148, 127. 

[73] Pachos J K et al 2008 Phys. Rev. Lett. 100, 156806. 

[74] Jarosik N et al 2010 Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Sky Maps, 

Systematic Errors, and Basic Results Preprint arXiv:1001.4744 

[75] Fulling S A 1973 Phys. Rev. D 7, 2850. 

[76] Davies P C W 1975 J. Phys. A 8, 609. 

[77] Unruh W G 1976 Phys. Rev. D 14, 870. 

[78] Takagi S 1986 Progress of Theoretical Physics Supplement 88, 1. 

[79] Hawking S W 1974 Nature 248, 30. 

[80] Singh D and Mobed N 2009 Class. Quantum Grav. 26, 185007. 

[81] Cortijo A and Vozmediano M A 2007 Nucl. Phys. B 763, 293. 

[82] Gibertini M et al 2010 Phys. Rev. B 81, 125437. 

[83] Verstraete F et al 2009 Phys. Rev. A 79, 032316. 

[84] Latorre J I 2010 Phys. Rev. A 81, 060309. 

[85] Susskind L 1977 Phys. Rev. D 16, 3031. 

[86] Rindler W 1966 Am. J. Phys. 34, 1 174. 

[87] Meyer J C et at 2007 Nature 446, 60. 

[88] Meyer J C et al 2007 Solid State Communications 143, 101. 



Dirac Equation For Cold Atoms In Artificial Curved Spacetimes 



20 



[89] Stolyarova E et al 2007 Proc. Natl. Acad. Sci. 29, 9209. 
[90] Ishigami M et al 2007 Nana Letters 7, 1643. 
[91] Bermudez A et al 2010 New J. Phys. 12, 033041. 

[92] Bakr W S et al 2009 A quantum gas microscope - detecting single atoms in a Hubbard regime optical lattice 
Preprinf larXiv:0908.0174l 

[93] Bakr W S et al 2010 Probing the Superfluid to Mott Insulator Transition at the Single Atom Level Preprint 
larXiv: 1006.07541 

[94] Sherson J F et al 2010 Single-Atom Resolved Fluorescence Imaging of an Atomic Mott Insulator Preprint 

larXiv:1006.3799l 
[95] Julsgaard B 2001 Nature 413, 400. 
[96] Eckert K et al 2007 Phys. Rev. Lett. 98, 100404. 
[97] Eckert K et al 2008 Nat. Phys. 4, 50. 
[98] Ernst P T et al 2009 Nat. Phys. 6, 56. 

[99] Wald R M 1984 General Relativity, First edition ed. (University Of Chicago Press). 
[100] IorioA2010 Weyl-Gauge Symmetry of Graphene Preprint \wXiv: 1007.50121 

[101] Mazza L et al 2010 Emerging Bosons with Three-Body Interactions from Spin-1 Atoms in Optical Lattices 

Preprint \arXiv: 1007.23441 
[102] Cortijo A and Vozmediano M A 2009 Nucl. Phys. B 807, 659, (Erratum). 
[103] de Juan F et al 2007 Phys. Rev. B 76, 165409. 
[104] Sitenko Y and Vlasii N 2007 Nucl. Phys. B 787, 241. 



